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1 .  INTRODUCTION 


This  document  comprises  the  final  report  on  the  Air  Force  Contract  no. 
F33615-84-C-3619  "Computer-Aided  Design  of  Robust  Decentralized  Controllers 
based  on  the  Stable  Factorization  approach.  For  the  Phase  I  (feasibility) 
effort  we  proposed  "to  make  a  start  on  the  design  of  robust  decentralized 
controllers  by  tackling  the  following  problems:" 

1)  numerical  computation  of  optimally  robust  controllers 

2)  decentralized  stabilization  of  an  n-channel  system 

3)  robust  decentralized  control  of  a  two-channel  system 

4)  optimal  decentralized  filtering  for  a  two-channel  system 

In  addition  to  the  Phase  I  proposed  tasks,  in  this  report  we  include  a 
new  subtask  ( 1 . ( b) )  on  simultaneous  stabilization.  Also  a  model  of  a 
flexible  structure  is  presented  for  which  two  types  of  "robust"  controllers 
are  computed. 

We  are  happy  to  report  that  substantial  progress  was  obtained  in  all 
of  the  four  original  tasks  as  well  as  the  added  substask  l.(a).  Our  preli¬ 
minary  computational  experience  with  the  design  of  the  above-mentioned 
controllers  has  been  very  good. 

As  a  result  of  successful  completion  of  Phase  I  effort,  a  firm  foun¬ 
dation  has  been  laid  for  the  Phase  II  undertaking. 


2.  TASK  1 


NUMERICAL  COMPUTATION  OF  OPTIMALLY  ROBUST  CONTROLLERS 

The  overall  goal  of  this  task  was  to  devise  numerical  techniques  for 
obtaining  optimally  robust  controllers.  The  specific  objective  of  this 
task  was  to  obtain  controllers  which  are  robust  against  the  class  of  plant 
perturbations  introduced  in  Vidyasagar  and  Kimura  (1984)  known  as  the 
"stable  factor  perturbations".  Recall  that  this  class  ,  denoted  as  S  (Ng  , 
Dq  ,  r  ),  consists  of  all  plants  P  of  the  form  P  (s)  =  N  (s)  [  D  (s)  ]”1 
for  which  | |  [  N'  -  N'g  D'  -  D'g  ]  j  I  <  |r(jw)  |  for  all  w,  where  r  is  a 
proper  and  stable  rational  function  and  prime  denotes  transposition.  Thus 
the  class  of  stable  factor  perturbations  consists  of  all  plants  whose 
stable  factors  are  'close'  to  the  stable  factors  of  the  nominal  plant  Pg. 
The  function  r(jw)  plays  the  role  of  an  "uncertainty  profile",  as  a  func¬ 
tion  of  frequency. 

The  advantage  of  stable  factor  perturbations  over  for  example 
multiplicative  or  additive  perturbations  is  that  the  former  does  not 
have  the  often  restrictive  requirement  on  the  perturbed  and  unperturbed 
plants  that  they  have  the  same  number  of  poles  in  the  unstable  region. 
Furthermore,  stable  factor  perturbation  results  do  not  require  that  the 
perturbed  plant  poles  or  zeros  on  the  jw-axis  be  fixed  in  their  location. 

In  order  to  illustrate  the  significance  of  these  more  relaxed  conditions,  a 
model  of  a  flexible  structure  will  be  discussed,  and  its  relevance  to 
stable  factor  uncertainty  will  be  described. 
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In  addition  to  the  tasks  spelled  out  in  the  proposed  Phase  I  effort, 
we  will  present  a  second  type  of  robustness  result.  (Furthermore  the 
robustness  results  will  be  applied  to  stabilize  a  flexible  structure.) 

This  second  type  of  robustness  relates  to  the  concept  of  simultaneous  sta¬ 
bilization  introduced  in  Vidyasagar  and  Viswanadham  (1982).  Recall  that 
the  simultaneous  stabilization  problem  is  as  follows:  Given  the  nominal 

plant  P0  and  a  number  of  contingent  plant  conditions  Pj  ,  -  ,  Pr  , 

determine  a  single  controller,  when  one  exists,  such  that  it  stabilizes  not 

only  the  nominal  plant  but  also  the  other  contingent  plants  P^  ,  - ,  Pr  . 

(For  details  see  Razavi,  Mehra,  and  Vidyasagar  (1984)  and  the  above  paper). 

Before  we  present  the  techniques  for  solving  the  above  two  types  of 
robustness  problems  we  describe  a  simple  model  we  shall  use  to  obtain  spe¬ 
cific  controller  designs.  Consider  a  4-disk  vibrating  system  indicated  in 
Figure  1. 
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The  disks  are  connected  by  torsion  springs.  Let  Kj.  denote  the  spring 
constant  for  the  spring  between  i-th  and  (i+l)-th  disk;  let  denote  the 


inertia  for  i-th  disk.  The  equations  of  motion  for  the  system  are: 
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where  Of  denotes  the  angle  of  deflection  of  the  i-th  disk  from  normal 
position.  Each  dot  indicates  differentiation  with  respect  to  time. 

The  right  hand  side  indicates  the  scaled  torques  applied  (in  this  case  to 
the  3-rd  disk). 


The  motivation  for  this  model  is  that  by  applying  actuator  inputs 
(torques)  uj  at  disk  i  one  wishes  to  control  the  outputs  0 j f  for  i  equal  to 
or  different  from  j.  The  transfer  function  from  each  input  u£  to  each  out¬ 
put  Oj  can  be  easily  calculated  from  the  dynamic  system  (S).  For  the  pur¬ 
pose  of  illustration  we  assume  that  the  spring  constants  Kj  are  all  equal 
to  unity.  The  nominal  plant  is  further  assumed  to  have  =  1  ,  i=l,  2,  3, 
4  (e.g.  equal  inertia).  Note  that  due  to  the  fact  that  there  is  no 


/. 


damping  in  the  system,  the  transfer  function  from  any  input  to  any  output 


will  have  all  of  its  poles  and  zeros  located  on  the  jw-axis. 

For  example  the  transfer  function  from  U3  to  ©4  will  be  of  the  form 


a  (s) 

S2  +  02 

=  _ 1 _ 

4  \ 

S2  .11,  (S2  +  ot2  ' 

i*l  i 

Where  3^  and  aj  are  real  positive  scalars,  corresponding  to  a  zero  at 
zj  =  and  poles  Pi  =  jail,  i=  1,  2,  3,  4.  The  rigid  body  modes 
correspond  to  the  poles  s2  *  0. 

The  preceding  model  will  be  used  in  the  application  of  two  robust 
controller  design  techniques  based  on  the  robustness  notions  that  we 
discussed  earlier:  (a)  robust  controllers  for  'discrete'  or  contingent 
plant  perturbations  (e.g.  fixed  controllers  that  simultaneously  stabilize  a 

nominal  plant  Pj  together  with  contingent  plants  P2  , _ ,  Pr.).  Here 

we  consider  the  one  input  -  two  output  case;  furthermore  let  r=2,  i.e.  one 
contingency,  (b)  optimally  robust  controllers  for  stable  factor  pertur¬ 
bations.  Sections  l.(a)  and  l.(b)  describe  the  computational  techniques 


together  with  numerical  results 


In  this  section,  we  present  an  example  of  the  simultaneous  stabiliza¬ 


tion  technique  introduced  in  Vidyasagar  and  Viswanadhara  (1982).  The 
technique  is  applied  to  a  four-disk  vibrating  system  described  earlier. 

The  computation  of  the  actual  controller  using  state-  space  versions  of  the 
stable  factorization  approach  is  carried  out  using  the  methods  from  the 
forthcoming  thesis  of  Minto  (1985). 

A  control  input  is  applied  to  the  second  disk  from  the  top,  and  the  two 
outputs  are  the  angular  positions  of  the  first  and  fourth  disks  respec¬ 
tively.  Two  different  situations  are  considered.  In  the  first,  all  four 
disks  are  identical.  In  the  second,  the  first  three  disks  are  identical 
and  the  last  disk  has  a  moment  of  inertia  which  is  one  quarter  of  those  of 
the  rest.  In  the  nominal  case  the  system  is  denoted  by  (s)  and  equals 
the  following:  (The  particular  pole  zero  locations  correspond  to  the 
experimental  set-up  described  in  Cannon  and  Rosenthal  (1984)) 

Case  1  (All  disks  identical) 
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Where  Che  symbol  (stuj)  denotes  (s+aj)(s-aj) 

Case  2  (First  3  disks  are  identical  and  moment  of  inertia  of  the  last 
disk  is  one  quarter  of  each  of  the  rest) 


P2(s) 


1 

a2  (s) 


f  b2i(s) 
|_b22(s) 


(1.3) 


where 


a2(s) 

=  s(s±j  29.34)  (s±j  21.48)  (s±j  11.938) 

( 1 .4a) 

b2i(s) 

=  b22<s) 

(1.4b) 

b22( s ) 

=  (s±j  25.95) 

(1.4c) 

Note  that  we  have  assumed  a  unity  gain  since  non-unity  gains  can  be  easily 
absorbed  into  the  controller.  The  objective  is  to  find  a  single  controller 
C  that  stabilizes  both  PI  and  P2.  In  the  present  instance  both  systems  are 
entirely  undamped,  l.e.  all  poles  and  zeros  are  purely  imaginary.  If  we 
stabilize  the  two  systems  vy  moving  the  poles  into  the  left  half-plane,  it 
is  possible  to  come  up  with  designs  that  yield  closed-loop  systems  which 
are  nominally  stable,  but  have  poles  with  very  low  damping.  To  rule  out 
this  possibility,  we  define  the  stability  region  to  be  the  region  (s:Re  s 
<-!}.  Thus,  in  order  for  a  system  to  be  stable  in  this  more  restrictive 
sense,  its  poles  must  have  real  parts  less  than  -1. 


Mathematically,  the  easiest  way  to  handle  this  requirement  is  to  trans¬ 
form  variables  by  replacing  s  by  s-1.  This  gives  two  modified  systems 


^  1  ^  ’CL  a~- »  -  •- 


•-  -  -t-  -t-  Q-  J 
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Q}(s)  and  (^(s)  whose  poles  and  zeros  all  lie  along  the  shifted  vertical 
axis  Re  s*l.  If  we  now  simultanteously  stabilize  Qj  and  Q2  by  a  controller 
Cmod(s)>  then  the  controller  C(s)“Cmo<j(s+l)  moves  the  poles  of  Pj  and  P2 
into  the  region  Re  s<-l. 

Accordingly,  define 


bl  1  (  s—  1 5“ 

1 

b2i(s-l) 

Q2  (s)= 

,  a2(s-l) 

bi2(s-l)_ 

i?22(8“l)  _ 

(1.5) 


where  all  functions  are  as  in  (1.2)  and  (1.4).  Because  the  systems  have 
only  one  input,  it  is  easy  to  construct  right-coprime  factorizations  for 
Ql  and  Q2  .  Define,  for  example 

f 1 ( s ) =( s+ 1 ) ( s+2 ) ( s+3 ) ( s+4 ) ( s+5 ) ( s+6 ) ( s+7 ) ,  (1.6a) 
f2(s)=(s+1.5)(s+2.5)(s+3.5)(s+4.5)(s+5.5)(s+6.5)(s+7.5)  (1.6b) 
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nil 


di  )is  an  r.c.p.  of  Qj  for  1*1,2.  It  is  easy  to  construct 


1— ni2  J 


state-space  realizations  for  the  matrices  in  (1.7).  Expand  the  various 


polynomials  as  follows: 


A 


Now  let  Af  be  the  companion  matrix 


-  4>i0  “  4>11 


“  ♦IS  “  <t>16 


r7x7. 


(1.9a) 


and  define 


0  ...  1 


Q  0  0  .  .  .  o]  gR7x1, 


(1.9b) 


C1  "  b110 


a10  “  a10  •  •  •  ®16  “  a16 


eR^x7 f 


(1.9c) 


Then  the  quadruple  (Ai,Bi,Ci,E^)  Is  a  state-space  realization  of  the 
"stacked"  3x1  system  [d^  n^j  nl2]’*  Here  E  has  been  used  in  place  of  the 
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>Svi>S>5>£4>5v3;«S-^>i»^ 


usual  output  feedthrough  symbol  D  to  avoid  confusion  with  'denominator' 
term  in  coprime  factors. 


Next,  consider  the  3x2  system 


dl  d2 


nll  n2 1  =:T 

,n12  n22_. 


A  state-space  realization  of  T(s)  can  now  be  readily  constructed,  as 


(1.10) 


follows: 


Ai  0 

It  * 

_0  A2_ 


eR14xl4, 


Bl  0 
_0  B2_ 


eR14x2 


Ct  -  [C1  c2]  e 1 4 f 


Et  *  [Ex  E2]  eR3x2 


(1.11a) 


(1.11b) 


Now  we  summarize  the  procedure  of  Vidyasagar  (1985),  Vidyasagar  et  al. 
(1985)  for  simultaneously  stabilizing  the  two  plants  Qj  and  Q2  .  Select 
units  u[  and  u2  in  S,  and  try  to  find  a  matrix  R(s)eSlx3  such  that 


R(s)  T( 8 )  -  (ui(s)  u2(s)]. 


(1.12) 


If  such  an  R  can  be  found,  then  R  can  be  partitioned  as  [dc  nci  nc2] , 


and  then 


1  fnci  (s)  nc-( s) ] 
C(s)  =  dc(s) 


(1.13) 
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Is  a  controller  that  solves  the  simultaneous  stabilization  problem.  The 
argument  in  Vidyasagar( 1985) ,  Vidyasagar  et  al.  (1985)  is  that,  since 
T  has  more  rows  than  columns,  generically  it  has  a  stable  left  inverse,  say 
Z.  Then  we  can  choose  the  units  m  and  U2  arbitrarily,  and  set  R=»[ui  U2]Z. 

However,  in  the  present  instance,  T  does  not  have  a  stable  left 
inverse,  since  the  matrix  T(  «)=Et  only  has  rank  1.  In  this  case,  we 
proceed  as  follows  (Note  that  the  procedure  is  from  Minto  (1985)): 

Let  t}(s),  t2(s)  denote  the  two  columns  of  T(s),  and  partition  B^E^as 

Bt  =  [bi  b2] ,  Et  »  [e!  e2].  (1.14) 

(Hereafter  we  drop  the  subscript  “t"  on  all  matrices.)  Then  it  is  easy  to 
see  that  the  quadruple  (A,bi,c,ei)  is  a  state-space  realization  of  t£. 

Also,  ei  =e2«  Now,  we  want  to  solve 

R  [q  t2]  =  [ut  u2].  (1.15) 

This  is  equivalent  to 

R  [ t  j  t2~t1]  -  [uj  u2-uj ] .  (1.16) 

Now  (A,b2~bj ,c,e2~e^)  is  a  state-space  realization  of  t2~tj.  Since 
e2“el=0,  t2(8)-ti(s)  is  strictly  proper.  Hence,  (1.16)  implies  that  U2 
-U1  is  also  strictly  proper.  In  other  words,  the  two  units  uj  and  U2  can 
not  be  chosen  independently,  but  must  satisfy  ^(“O^u^®).  Now  (1.16)  is 
equivalent  to 


R  [tt  (s  +  y)  (t2  ~t!  )]  *  [ u j  (s  +  y)  (u2  -  uj  )], 


(1.17) 


where  y  Is  any  positive  number.  Since  (A,b2  -  bj,  c,0)is  a  state-space 
realization  of  t2-tj,  one  can  show  without  much  difficulty  that  a  state- 
space  realization  of  (s+y )(t2(s)-tj(s))  is  given  by  (A,(A+Yl)(b2-bj) ,c, 
c(b2-bj)).  If  the  vector  c(b2-bj)  is  linearly  independent  of  e^,  then  the 
matrix  [tj  (s+y )( t2-tj ) ]  has  a  left  inverse,  and  (1.17)  can  be  solved  for 
R.  Otherwise  the  procedure  has  to  be  repeated.  The  algorithm  can  be 
stated  as  follows: 


Step  1:  Set  i=*0 

Step  2:  Is  e2  linearly  independent  of  e^  ?  If  so,  go  to  step  5; 
if  not,  go  to  step  3 

Step  3:  Since  e2  is  linearly  dependent  on  ej,  there  exists  a  constant 

v^  such  that  e2  =  v^  ej  .  Replace  b2  by  b2  -  v^  bj  ,  e2  by  e2 
-  vj  ej  =  0. 

Step  4:  Replace  b2  by  (A+'yI)  b2  ,  e2  by  c  b2  .  Increment  i  and  go  to 
step  2. 

Step5:  Stop 

At  this  stage,  one  gets  an  equation  of  the  following  form 


R  (s)  [ti  (s)  t2  (s)  ]  =  [u2  (s)  u2  (s)  ], 


(1.18) 


where,  letting  k  denote  the  total  number  of  iterations,  we  have 


U2  (s)  =  U2  (s)  -  Z  v*  (s  +  y  )*  (8) 

i=0 


(1.19) 


Also,  the  termination  criterion  ensures  that  the  matrix  [tj  (s)  t2  (s)] 
has  a  stable  left  inverse.  Hence,  once  suitable  units  uj  and  U2  are 
chosen,  the  matrix  R  (s)  can  be  found  as 


R  (s)  =  [uj  (s)  U2  (s)  ]  Tl  (s),  (1.20) 


where  T1  (s)  is  a  stable  left  inverse  of  [tj  (s)  t2  (s)  ]. 


To  complete  the  discussion,  it  only  remains  to  show  how  to  choose  the  units 
ui  and  U2  .  Now  any  function  f  in  s  can  be  expanded  in  a  power  series  of 
the  form 


f  (s)  =  Z  f±  (s  +  y  )  _i  •  (1.21) 

i=0 


Now,  the  power  series  of  the  units  uj  and  U2  must  satisfy  certain 
relationships  to  ensure  that  various  quantities  are  strictly  proper.  For 
instance,  if  we  simply  choose  uj  (s)  =1,  then  the  first  k  terms  in  the 


power  series  of  U2  (s)  must  be  vq  +  V2  (s  +  y)  +•••+  v^-i  (s  +  y)“^+^» 
where  the  constants  vj  are  generated  in  step  3  of  the  algorithm.  Thus  any 

unit  U2  can  be  chosen,  so  that  the  above  condition  is  satisfied. 
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In  the  problem  at  hand,  we  chose  y 


8  ,  and  It  turned  out  that  the 


algorithm  had  to  be  run  f ive  times,  and  the  constants  were 

v0  “  1»  vi  “  -  3.5,  v2  =  461.58,  V3  =  6769.9,  V4  -  -75847  (1.22) 

If  we  perform  the  bilinear  transformation  z  =  (s-8)  /  (s+8),  which  maps  s=<» 

into  z=l,  we  conclude  that  the  power  series  (in  terras  of  z)  for  u2  must 

look  like 

u2  =  3.0798  +  3.7164  z  -  0.18256  z2  -  2.9765  z 3  -  1.1573  z4  +...  (1.23) 

Fortunately,  it  turns  out  that  if  we  add  just  one  term,  namely  10z5,  then 
u2  is  a  unit,  i.e.  its  zeros  are  in  the  open  left  half  of  the  s-^plane. 

The  procedure  leads  to  a  nineteenth  order  controller,  whose  state-space 
realization  is  given  in  Appendix  B.  The  twenty-six  closed-loop  poles  with 
plants  Pi  and  P2  are  given  below: 

With  Plant  P^  and  controller 

-1 

-1.517 

-1.5 

-1.333 

-1.25 

-1.20 

-1.167 


-  14  - 


1  ><%  >'*  »  •  V"»‘  /•  '•  "«  , '«  ”• 


t 


-1.143 


h 


V. 


-1.118 

-1.126 

-1 . 155+jO. 1236 
-1 .006jJ0.077 
-1 .002+J0.0048 
-1.002+^0.0044 
-1 .01 1+J0.004 
-1.0124J0.006 
-1.013+J0.006 
-1.013+J0.003 


With  Plant  P2  and  controller 


-2.481 

-1.667 

-1.317 

-1.4 

-1.286 

-1.222 

-1.043 

-1.095 

-1.018 

-1.017 

-1.015 

-1.013 
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‘  -1.016hJ0.124 

1  -1.006+J0.077 

j  -l.002jJ0.048 

-l.002jJ0.044 
-l.003jJ0.039 

,  -1.137+j0.037 

* 

[  -1.105HJ0.025 

i 

It  can  be  noted  that  all  closed-loop  poles  have  damping  ratios  very 
!  close  to  1. 


In  summary,  in  an  actual  controller  design  problem  one  would  apply 
order  reduction  techniques  to  replace  the  nineteenth  order  controller  by 
something  smaller.  One  would  also  scale  the  A,B,C,E  matrices  of  the 
controller  in  such  a  way  that  the  numbers  do  not  have  a  very  large  dynamic 
range . 

Subtask.  1(b):  Numerical  Computation  of  Optimally  Robust  Controller 

Recall  that  in  this  task  the  class  of  perturbations  that  we  are 
interested  in  are  those  of  'stable  factor';  this  class  was  denoted  by  S(Nq, 
Dq,  r).  The  following  necessary  and  sufficient  result  was  derived  by 
Vidyasagar  and  Kimura  (1984).  A  Controller  C  that  stabilizes  the  nominal 
plant  Pq  NqDo-!  also  stabilizes  all  plants  in  the  class  S(Nq,  Dq.t)  if 
and  only  if  1 [X  Y]rl<l,  where  C  =  -Y~^X  and  X  and  Y  are  the  solutions  to  the 
Bezout  identity.  Stated  more  generally,  given  the  class  S(Nq,  Dq,  r)  of 
stable  factor  perturbations,  there  exists  a  controller  which  robustly  sta- 


blitzes  all  plants  in  S(N0,  Dq,  r)  If  and  only  if  I [Y  X] '  +  f-N  5] 'RK1 
for  some  stable  rational  matrix  R,  where  (R,  B)  denote  the  right  coprime 
factorization  of  P.  Given  an  uncertainty  profile  reS,  consider  Xr  for  all 
real  values  of  the  parameter  X.  An  optimally  robust  controller  is  that 
which  stabilizes  the  plants  in  the  set  S(Nq,  Dq,  Xr)  with  the  largest  value 
of  X. 

With  this  review,  our  aim  in  this  subtask  is  to  compute  the  following 
minimization  problem: 

Min  It  F-GRH«  (1.24) 

ReH°° 

nxl  . 

where  Ge  Ra  ,  He  Ro0lcxm  with  n>l ,  k<m 

(other  cases  are  simpler  and  follow  from  this  problem.)  To  see  that  the 
above  minimization  solves  the  robust  controller  problem  described  above 
simply  make  the  associations: 

F  =:  [Y  X]'  (1.25) 

G  -:  [ -N  D] ’  (1.26) 

H  *  I  (1.27) 

The  robust  controller  is  obtained  from  the  optimal  R,  via  Youla's 
parameterization. 

We  will  call  the  above  minimization  problem  the  "tall"  plant  minimiza¬ 
tion  problem  as  the  free  variable  R  appears  on  the  right  of  G,  and  G  is 
"tall".  The  proposed  procedure  for  solving  (1.24)  turned  out  to  be  dif¬ 
ficult  to  implement  -  as  the  proposed  minimax  problem  became  intractable. 
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Instead  we  present  the  following  approach.  (Throughout  we  assume  that  G 
and  H  have  full  column  or  row  rank  as  appropriate;  otherwise  the  procedure 
can  be  slightly  modified.)  The  procedure  will  be  summarized  without  proof 
as  an  algorithm.  The  main  idea  is  that  by  a  series  of  manipulations  we 
will  transform  the  tall  plant  problem  (1.24)  to  a  square  plant  problem 
(1.28)  which  then  can  be  readily  solved. 

Step  1 : 


Compute  G  =  Gq»  (Gq  square) 


H  =  Hq'  Hji’,  (Hq  square) 


where  the  subscripts  I  and  0  stand  for  inner  and  outer,  respectively 


Let 


Q  -  Gq  R  Hq 


Step  2:  Compute  complementary  inner  for  Gj 


i)  Perform  row  permutations  if  necessary  to  get 


Gi 


M 


with  D  *  0 


ii)  Let  P(Z)  =  -[D'(Z_1)]_1  N(Z_1) 
Compute  r.c.f.  of  P 

p(z)  =  u(z)  vkz) 


mx(n-m) 


(n-m)x(n-m) 


Let  W  =  |_  vj  t  Rnx(n-m) 


Compute  inner  outer  factorization  of  W: 

w  =  w*  w0 


,,  (n-m)x(n-m)  nx(n-m) 

W0  c  IL  ,  Wt  o  R 

is  the  complementary  inner  for  G^: 


*  =  fGi:Wi] 

Where  »  c.  R  nxn  is  inner 


Step  3: 


Compute 

A 

B 

<p*F  =: 

_  C 

D  _ 

(Partition  v*F  such  that  Q  and  A  are  additlvely  compatible.  Note  that 
«*=:♦' (-s) . ) 


Step  4: 


Let  |i  *  max  |J  C  D 


Step  5: 

Choose  y  >  m 

Let  K  =*  [C  D] 

Step  6: 

Compute  L  =  (y2  I  -  K*  K)-1/2 


(Notation:  If  A  =  V^v,  then  V=:A^2) 


M  =  (I-BLL*B*)-1/2 

KB  Llloo  >  1  increase  y  and  go  to  Step  6,  otherwise  go 
to  Step  7 

Compute  J*  =  Min  IIMAL  -  SR* 

S  e  M(Hoo) 

i)  Let  A2  =  MAL;  A2  e  M(L«,) 

ii)  Compute  I.c.f.  of  A2, 

A2  =  (T(z))-1  U(z) 

iii)  Compute  inner  outer  factorization  of  T: 

T  =  Ti  Tq;  where  Tq  is  a  unit 

Then 

A2  =  (T^z))-!  V(z) 

Where 

V(z)  =  Tq1  U 


Then 


b  =  |  T±  1 

IIMAL  -  SH«  =  »T?d^  V  -  b  Sil„ 
adj 


where 


x 


indicates  adjoint  of  x. 


Thus  the  minimization  problem  (1.24)  becomes 

J*  =  Min  HTfdjV  -  b  SI.  (1.28) 

S  e  M(Hoo) 

The  right  hand  side  of  (1.28)  is  a  'square'  plant  minimization 
problem  which  can  be  solved  using  methods  of  Section  7.4  of 
Vidyasagar  (1985) 


Step  9: 

If  J*<1  decrease  y  and  go  to  Step  6 
If  J*>1  Increase  y  and  go  to  Step  6 
If  J*=l  go  to  Step  10 

Step  10: 

Set  Q  =  M"1  S  L_1 

then,  R*  =  Go"1  Q  H“l0 

is  the  optimum  R  in  (1.24). 

We  conclude  this  section  with  a  controller  design  which  achieves 
stable  factor  robustness  against  pole/zero  movements  on  the  jw-axis  of  the 
four  disk  vibrating  (non-minimum  phase)  system.  The  controller  order  turns 
out  to  be  smaller  (than  that  of  the  previous  controller)  but  some  of  the 
closed  loop  poles  are  not  as  well  damped. 


V*.  V  V  V', 


In  the  interest  of  brevity,  only  the  final  results  are  summarized 
here.  Nominal  plant  transfer  function  from  3r<*  to  4^  disk: 


P0(s)  =  k 


s2  +  195 

s2(s2+156)(s2+361)(s2+576) 


One  such  robustly  stabilizing  controller  c(s)  for  Pg(s)  is: 

(s-132.44)(s+27.14)(s+90.26)(s+3.6)(s+0.99)(s2-17.5  s  +  398.5) 


c(s)  = 


( s+190.35) ( s+37. 70)(s+0.489±30.935j )(s+l 1 .61 ±20.89 j )(s+8.86±5.307j ) 


The  value  of  the  gain  k  is  827,000  (Cannon  and  Rosenthal  (1984)).  The 
closed  loop  nominal  system  poles  are 


-1.9036  +  0.00001 
-0.0063  +  0.30961 
-0.0063  -  0.30961 
-0.0036  +  0.24271 
-0.1183  +  0.21201 
-0.0036  -  0.24271 
-0.1183  -  0.21201 
-0.0040  +  0.17721 
-0.0040  -  0.17721 
-0.0017  +  0.12351 
-0.1001  +  0.06031 
-0.0017  -  0.12351 
-0.1001  -  0.06031 
-0.0370  +  0.00001 
-0.0181  +  0.00001 
-0.0073  +  0.00001 

The  variations  in  inertia  cause  the  zero  and  the  poles  to  vary  from  their 
nominal  values.  Figure  2  indicates  the  range  for  the  zero  and  one  pole 
variation.  The  nominal  value  for  the  zero  Is  13.96  and  for  the  pole  is 
12.48.  The  acceptable  range  for  pole  zero  variations  is  shown  by  the 
darkened  lines.  Other  pole  variations  of  about  the  same  percentage 
variation  range  turn  out  to  be  stabilized  by  the  fixed  controller. 
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3.  TASK  2 


1 


SR© 


DECENTRALIZED  STABILIZATION  OF  AN  n-CHANNEL  SYSTEM 


Recall  that  task  2  consisted  of  finding  a  parameterization  of  all 
decentralized  controllers  that  stabilize  a  given  plant. 

Suppose  a  plant  P  has  dimensions  Jlxm,  so  that  it  has  m  inputs  and  1 
outputs.  Partition  the  Inputs  and  outputs  into  k  disjoint  subsets  each, 
where  k  is  the  number  of  channels.  Without  loss  of  generality,  renumber 
the  inputs  and  outputs  in  such  a  way  that  the  first  channel  consists  of  the 
first  mj  inputs  and  outputs,  the  second  channel  of  the  next  m2  inputs  and 


outputs,  and  so  on.  Then  we  have  integers  mj 


1%  ,  and  9^  j  , . . . , 


k  k 

fck  such  that  £  m^  =  m  and  £  =  l.  This  partitioning  of  inputs  and 

i=l  1=1 

outputs  is  called  an  information  structure.  With  these  definitions,  a 
controller  C  is  said  to  be  decentralized  if  it  is  block-diagonal  of  the 


Cl 

c=  c2 


(2.1) 


where  cj  has  dimensions  ra^xi^.  The  objectives  of  this  task  are  (i)  to 
derive  conditions  under  which  a  given  plant  can  be  stabilized  by  a 
decentralized  controller,  and  to  find  an  expression  for  'll  controllers 
that  stabilize  it. 


.  -  «  •  - 
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We  now  summarize  some  results  from  Vidyasagar  (1985)  concerning 
stabilization.  Given  a  plant  p  and  a  controller  c,  we  form  the  composite 


transfer  matrix 


H  (P,  O- 


Tl+PC)-1 

£(I+PC)_1 


-P(I+CP)-1 

(I+CP)_1_ 


(2.2) 


and  say  that  C  stabilizes  P  if  the  above  transfer  matrix  is  stable.  Let 
S  denote  the  set  of  scalar  stable  transfer  functions.  Thus  in  the  con- 
tineuos-time  case  S  consists  of  proper  rational  functions  of  the  Laplacian 
variable  s  whose  poles  all  lie  in  the  open  left  half-plane.  In  the 
discrete-time  case,  S  consists  of  proper  rational  functions  of  the  unit 
advance  operator  z  whose  poles  all  lie  inside  the  closed  unit  disk. 
Finally,  let  M  (S)  denote  the  set  of  matrices,  of  whatever  order,  whose 
elements  all  lie  in  S.  Then,  given  any  P,  one  can  factor  it  in  the  form 
Np  Dp1-  Dp1  Np,  where  Np,  Dp  "Np,  Dp  all  belong  to  M  (S),  and  moreover 
satisfy  the  so-called  Bezout  identities 


Xp  Np  +  Yp  Dp  -  I,  Np  Xp  +  Dp  Yp  -I,  (2.3) 

where  Xp,  Yp,  Xp,  Yp,  also  belong  to  M(S).  We  refer  to  the  pair  (Np,  Dp) 

as  a  right-coprlme  factorization  (r.c.f.)  of  P,  and  to  the  pair  (Dp,  Np) 
as  a  left-coprlme  factorization  (l.c.f.)  of  P.  Similar  remarks  and  symbols 
apply  to  C. 


A  function  f  in  S  is  called  a  unit  if  its  reciprocal  1/f  also  belongs 
to  S.  A  square  matrix  U  in  M  (S)  is  called  unimodular  if  its  inverse  U-1 
also  belongs  to  M  (s).  It  is  easy  to  show  that  U  is  unimodular  if  and  only 
if  its  determinant  (denoted  by  |u|)  is  a  unit. 


The  following  theorem  Is  well-known,  and  forms  the  basis  of  this 
section;  see  (Vidyasagar  1985). 

Theorem  2.1.  Given  P,  C,  let  (Np,Dp),  (NC>DC)  be  r.c.f.'s  of  P  and  C 
respectively,  and  let  (Dp,Np),  (DC,NC)  be  l.c.f.'s  of  P  and  C  respectively. 
Then  the  following  sattements  are  equivalent: 

(i)  C  stabilizes  P. 

(ii)  The  matrix  A  =  Dc  Dp  +  Nc  Np  is  unimodular. 

(iii)  The  matrix  A  =  Dp  Dc  +  Np  Nc  is  unimodular. 

Now  we  are  in  a  position  to  present  some  new  results.  A  key  concept 
in  decentralized  stabilization  is  that  of  fixed  modes  (Davison  and  Wang 
(1973),  Vidyasagar  and  Viswanadham  (1983).  In  the  present  report,  we 
extend  this  concept  by  defining  a  decentralized  fixed  determinant.  Given 
a  plant  p  and  a  controller  c,  it  can  be  shown  that  the  quantities  |a|  and 
|Z|  (where  A,  A  are  defined  in  theorem  2.1)  differ  only  by  a  unit,  and  that 
each  is  a  characteristic  determinant  of  the  closed-loop  transfer  matrix 
H(p,c);  see  (Vidyasagar  1985.).  The  significance  of  this  is  that 
the  unstable  poles  of  H(p,c)  are  precisely  the  RHP  zeros  of  its 
characteristis  determinant.  Now  the  (decentralized)  fixed  determinant  <|>  is 
defined  as  follows:  Let  <j>  =  <f>(p,c)  denote  the  quantity  |a|,  and  let  R 
denote  the  set  of  constant  decentralized  feedback  gain  matrices.  Then  f, 
the  fixed  determinant,  is  given  by 

f  =  g.c.d.  4>  (p,K) 

(2.3) 

K  c  R 

In  other  words,  f  is  the  greatest  common  divisor  of  all  characteristric 
determinants  of  all  closed-loop  systems  obtained  by  applying  a  block- 


diagonal  constant  output  feedback  to  the  given  plant  P.  Obviously  the  RHP 
zeros  of  P  are  precisely  the  RHP  poles  of  P  that  remain  fixed  under  all 
such  feedback  gains.  The  next  two  theorems  give  some  new  and  powerful 
results  characterizing  f.  (All  proofs  are  given  in  Appendix  A.) 


Theorem  2.2.  Let  (N,D)  be  an  r.c.f.  of  P,  (D,  N)  an  l.c.f.  of  P.  Then 
each  of  the  following  constructions  gives  an  explicit  formula  for  f:  (i) 
partition  N  and  D  by  rows,  as  follows:  consists  of  the  first  m^  rows, 

of  D,  D2  of  the  next  m2  rows,  and  so  on  until  Dk  consists  of  the  last  1% 
rows .  Thus 


Dk 


where  D1  has  mj  rows.  Similarly  partition  N  as 


r  ni 

i  N2 


L  Nk  J 

where  N*  has  lj  rows. 


Now  define 


"  D1  “• 

N1  J 


and  note  that  F  has  dimensions  (m+H)  xm.  Then  f  is  the  g.c.d.  of  all  raxm 
minors  of  F  consisting  of  exactly  mj[  rows  from  F^,  for  each  i. 


(il)  partition  D,  N  as 


where  has  columns,  Nj  has  columns.  Define 

F±  =  [D*  Nil,  F  =  [Fj  ?2  ...Fk],  (2.8 

and  note  that  F  has  dimensions  £x[£+ra].  Then  f  Is  the  g.c.d.  of  all  SLxl 
minors  of  F  consisting  of  exactly  columns  from  F^. 

Theorem  2.2  gives  a  very  simple  and  direct  characterization  of  fixed 
modes,  or  more  generally  fixed  determinents.  Theorem  2.3  below  shows  that 
an  RHP  pole  of  P  that  can  not  be  moved  by  constant  (decentralized)  output 
feedback  also  remains  fixed  under  dynamic  (decentralized)  output  compen¬ 
sation.  Strictly  speaking  this  result  is  not  new;  however,  the  proof  given 
here  is  very  transparent. 

Theorem  2.3.  Let  DC  denote  the  set  of  all  decentralized  controllers 

and  define 

g  =  g.c.d.  4>(P,C)  (2.9) 

C  e  DC 
Then  g  =  f. 

Next,  we  present  a  decomposition  principle  that  is  useful  in  finding 
all  decentralized  stabilizing  controllers.  This  result  should  be  compared 
with  (Zames  1981),  (Callier  and  Desoer  1982)  and  (Vidyasagar  1985). 

Theorem  2. 4. Given  a  plant  P,  let  DS  (P)  denote  the  set  of  all 
decentralized  stabilizing  controllers  that  stabilize  P.  Let  C  be  a 
specific  decentralized  stabilizing  controller  for  P,  and  define  Pj*P 


( I+CP)-1 .  Then  C+DS  (P^)  is  a  subset  of  DS  (P),  and  the  two  sets  are 
equal  If  and  only  if  C  is  stable. 

The  theorem  should  be  interpreted  as  follows:  Given  a  plant  P,  suppose 
one  finds  first  a  particular  C  that  stabilizes  P.  Then,  whenever  Cj 
stabilizes  the  resulting  plant  Pj  =»  P(I+CP)-^,  the  controller  C+C^ 
stabilizes  the  original  plant  P.  If  C  is  stable,  then  every  stabilizing 
controller  for  P  can  be  expressed  as  the  sum  of  C  and  a  stabilizing 
controller  for  Pj.  If  C  is  unstable,  then  there  exist  stabilizing 
controllers  for  P  that  can  not  be  decomposed  in  this  fashion.  The  point 
of  the  theorem  is  that  the  set  DS  (Pj)  is  in  general  easier  to  determine 
than  the  set  DS  (P),  since  the  first  plant  is  stable.  Thus,  by  using 
Theorem  2.4,  one  can  reduce  the  problem  of  finding  all  decentralized 
stabilizing  controllers  for  a  general  plant  to  one  of  finding  all 
decentralized  stabilizing  controllers  for  a  stable  plant,  under  certain 
conditions.  Specifically,  if  the  given  plant  can  be  stabilized  by  a  stable 
decentralized  controller,  such  a  reduction  is  always  possible. 

Youla  et  al.  (1974)  call  a  plant  strongly  stabilizable  If  it  can  be 
stabilized  by  a  stable  controller.  Theorem  2.4  brings  out  the  importance 
of  the  concept  of  decentralized  strong  stabilizability ,  which  is  the  abi¬ 
lity  to  be  stabilized  by  a  decentralized  stable  controller.  While 
necessary  and  sufficient  conditions  are  available  for  strong  stabilizabi¬ 
lity,  so  far  only  necessary  conditions  are  available  for  decentralized 
strong  stabilizability. 

At  this  point,  it  becomes  necessary  to  choose  between  stating  all 
results  in  their  most  general  form  or  in  their  clearest  form,  we  opt 
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to  state  all  subsequent  results  in  their  clearest  form  by  restricting 


u 


to  the  case  where  the  plant  P  is  square  of  dimension  mxra,  and  there  are 
m  channels.  Thus  a  decentralized  controller  in  this  context  is  a  diagonal 
controller.  It  must  be  emphasized  that  all  of  the  results  below  carry  to  a 
general  plant  with  a  general  information  structure,  but  at  the  expense  of 
more  cumbersome  notation. 


Let  (N,D)  be  an  r.c.p.  of  the  plant  P.  Then  both  N  and  D  are  square 
and  have  dimension  raxm,  and  the  matrices  D*,  N*  in  (2.4)  -(2.5)  are  just 
the  rows  of  D  and  of  N.  Hence  we  choose  to  denote  them  by  lower  case 
letters  d^,  n*.  The  matrix  F  in  (2.6)  is  2mxm  and  consits  of  an 
interlearning  of  the  rows  of  D  and  N.  The  minors  of  F  whose  g.c.d,  equals 
the  fixed  determinant  f  are  those  of  the  form 


r  vi 
I  v2 


(2.10) 


where  v*  equals  either  d1  or  n*.  Clearly  there  are  2m  of  these.  An  easy  way 
to  index  these  minors  is  to  use  binary  notation.  Let  a  e  {0,l}m,  i.  e.  a 
is  a  binary  number  with  m  bits.  Let  v^=  d*  if  =  0  and  n*  if  otf  =  1. 

Then,  as  o  varies  over  { 0 , 1 } m ,  the  corresponding  minor  va  defined  in  (2.10) 
traces  out  all  the  minors  whose  g.c.d.  equals  the  fixed  determinant  f. 

Theorem  2.5.  Let  ()  denote  the  binary  zero  with  m  bits,  In  order  for 
P  to  be  stabilizable  by  a  stable  diagonal  controller,  a  necessary  condition 
is  the  following.  Define 


>  V. 


(2.11) 


-»»rJ 


a  “  v£  =  |D| ,  b  =  g.c.d.  va 
a  *  0 

Then  a  and  b  must  satisfy  the  parity  interlacing  property;  that  is,  between 
every  pair  of  real  RHP  zeros  of  b,  there  must  be  an  even  number  of  zeros  of 


a,  or  equivalently,  a  must  have  the  same  sign  at  all  real  RHP  zeros  of  b. 


The  condition  of  Theorem  2.5  should  be  compared  to  the  condition  for 
strong  stabilizability  by  a  (not  necessarily  decentralized)  controller,  as 
given  in  Youla  et  al  (1974),  Vidyasagar  and  Viswanadham  (1982):  Let  n 
denote  the  g.c.d.  of  all  elements  of  N;  then  a  and  n  satisfy  the  parity 


interlacing  property.  The  latter  condition  is  necessary  and  sufficient, 
while  the  former  (in  Theorem  2.5)  is  only  necessary. 


Next,  we  turn  to  the  problem  of  finding  all  decentralized  stabilizing 
controllers.  If  the  plant  is  strongly  decentrally  stabilizable ,  then  the 
results  below  can  be  simplified  using  Theorem  2.4.  The  following  lemma 
is  very  useful  for  this  purpose. 


Lemma  2.6.  Let  gij  l<i<£,  0<j<m  be  given  elements  of  S,  and  suppose 


that  for  some  integer  j>l  the  matrix 


f  810 


(2.12) 


g*jJ 


has  rank  2.  Define 

a  =  g.c.d.  { gi j }  (2.13) 

1  <  i  <  0  <  j  <  m 

under  these  conditions,  for  almost  all  functions  a j ,  l<j<m  belong  to  S, 


we  have 


m 

g.c.d.  {gl0  +  E  gti  ai }  =a  (2.14 

j-1 

l<i<l 

Now  we  can  tackle  the  problem  of  finding  all  decentralized  stabilized 
controllers  using  the  results  of  Corfmat  and  Morse  (1976),  together  with 
Lemma  2.6.  We  first  give  a  detailed  treatment  of  a  two-input,  two-outputs 
systems  to  bring  out  the  basic  ideas  clearly,  and  then  move  on  to  the 
general  case  of  an  m-input,  m-output  system. 


consider  a  system 


PI  2 

P22 _ 


which  we  wish  to  stabilize  by  a  diagonal  controller 


(2.15) 


C1 

0 


0 


(2.16) 


Let  (D,FJ)  be  an  l.c.f.  of  P,  and  partition  the  matrix  [0  N]  as 
[D  FI]  =  [dt  d'2  n j  n2] 


(2.17) 


Since  C  is  diagonal,  it  ras  an  r.c.p.  of  the  form 


E 


just  let  (8^,04)  be  a  coprime  factorization  of  c^,  for  i=l,2.  In  order  for 
C  to  stabilize  P,  by  theorem  2.1  a  necessary  and  sufficient  condition  is 
that  the  quantity 


u=  |DDc  +  Nnc| 

be  a  unit  of  S.  However,  from  (2.17)  and  (2.18), 


(2.19) 


DDC  +  Nnc  =  [0  FJ] 


"Dcl  ‘  [ 

-»c- 


04  dj  +  8i  nj  a2  d2  +  $2  n  2 


+  8?  n  2 -i  (2.; 


Now  for  convenience  we  introduce  the  "wedge  product"  notation. 

Given  two  2x1  vectors  vj  and  V£,  the  scalar  vy  A  V2  equals  the  determinant 
of  the  2x2  matrix  [v!  V2] .  It  is  easy  to  see  that  wedge  product  is  a 
multilinear  function.  Now 


u=  (04  dj  +  8^  nj)  A  (<*2  d2  +  82  n2^ 

=  04  a.2  (diAd2)  +  04  82  (djA^)  +  81  04  (nlAd2^  +  Pi  P2  (niA°2) 


(2.21) 


clearly,  unless 


g.c.d.  { d 1 Ad2 »  diAn2,  n^Ad2»  njAn2)  =1, 


(2.22) 


u  can  never  be  a  unit.  It  can  be  routinely  verified  that  the  above  g.c.d. 


is  precisely  the  fixed  determinant  defined  in  Theorem  2.2 


Suppose  (2.22)  is  true.  Then  the  problem  of  choosing  a  decentralized 
stabilizing  controller  becomes  one  of  selecting  04,  <*2 ,  Si,  02,  such  that 
u  is  a  unit.  For  this  purpose,  note  that 

u=  ai  [<*2  (d iAd2)  +  02  (djA^)]  +  01  Ea2 

+  02  ( n i A n2 ) J  (2.23) 

Hence,  if 

g.c.d.j  a2  (d[Ad2)  +  02  (djAn2),  02  (nl^n2^  +  02  ( n l ^n2 ) }  =1 «  (2.24) 

then  it  is  possible  to  choose  04,  0i  such  that  u  is  a  unit;  the  converse 
is  also  true.  Thus  the  question  becomes:  When  is  (2.24)  true  ?  This  is 
where  Lemma  2.6  is  useful.  It  states  that  if  (2.22)  holds,  then  (2.24) 
also  holds  for  almost  all  02*  02»  provided 


d  L  Ad2 
ni  Ad2 


d!An2 

ni  An2 


g*0 


However,  routine  calculations  show  that 

g  =  |  D  j  P12  P21  (2.26) 

Hence  the  rank  condition  holds  provided  P12  +  0  and  P21  +  0.  This  is  the 
strong  connectedness  condition  of  comfort  and  Morse  (1976). 

Thus  the  strategy  for  finding  all  decentralized  stabilizing  controllers  is 
as  follows:  Pick  02,  02  at  random.  Then  for  almost  all  choices  of  02, 02 
the  condition  (2.24)  will  hold.  Next,  select  x,y  such  that 

x  f '12  ( fl  1  Ad 2 )  +  02  (djAn2)j  +  y  [02  (njAd2)  +  02  (nlAn2)]  =  1 

(2.27) 
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Then  the  set  of  all  04,  8j  such  that  u=l  (and  this  is  the  only  case  that 
needs  to  be  considered)  is  given  by 

04  =  x+r  [<*2  (niAd2)  +  82  (niAn2)]> 

81  =  y-r  [a2  (djAd2)  +  82  (djAnj)],  (2.28) 

as  r  varies  over  the  set  S.  The  formula  (2.28)  can  be  made  cleaner  by 
noting  that  P*B_i  &,  and  using  the  relationships  between  minors  of  p  and 
minors  of  [  8  R  ]  (see  Vidyasagar  1985,  appendix  B)).  This  gives 

04  =  x  +  r  |B|  [a2  Pll  +  02  pj2  ],  (2.29) 

8l  =  y  -  r  |D | [o2  +■  82  P22]. 

The  above  stated  discussion  becomes  a  little  simpler  if  the  plant  P  is 
already  stable.  As  shown  by  Theorem  2.4,  it  is  possible  under  certain  cir¬ 
cumstances  to  restrict  one’s  own  attention  to  this  case.  To  proceed 
further,  we  need  the  following  result: 

Theorem  2.7  (Youla  et  all  1976,  Zames  1981)  given  a  stable  plant  P, 
the  set  of  all  controllers  that  stabilize  P  is  given  by 

s(P)  .  ((I-RP)-l  R:  ReM<8>>  <2-30 

Using  this  theorem,  we  can  characterize  all  decentralized  stabilizing 
controllers . 

Theorem  2.8  Suppose  P  is  a  given  2x2  stable  plant.  Then  the  set  of 
all  diagonal  C  that  stabilize  P  is  given  by 

DS  (P)  =  {(I-RP)-1R:  ReM(s)  and  rl2  +  Pi 2  |  R  |  =  0,  r21  +  p21  |  R  |  =  0} 

(2.31) 


Thus  any  stable  R  substituted  into  the  expression  (I-RP)_j  R  yields  a 


controller  that  stabilizes  P.  However,  in  order  for  the  resulting 
controller  to  be  diagonal,  R  must  satisfy  the  two  additional  conditions 
r12  +  p12  |  R|  =  0,  r2\  +  P2i  |  R  |  =  0.  These  conditions  can  be  used  as 
follows:  pick  any  function  f  e  S,  and  set  r^2  =  “Pl2  r21  = 

-p2i  f,  and  select  r^,  r^  arbitrarily  except  for  the  constraint  r^ 

r22  5,1  f  0-P12  P21  f)* 


4.  TASK  3 


ROBUST  DECENTRALIZED  CONTROL  OF  A  TWO-CHANNEL  SYSTEM 

An  important  practical  problem  is  the  following:  Given  a  nominal 
plant  description  Pg,  together  with  a  "band  of  uncertainty"  containing  Pq, 
when  does  there  exist  a  controller  that  stabilizes  not  only  Pq  but  also  all 
plants  within  the  band  of  uncertainty?  This  problem  can  be  divided  into 
two  parts:  (i)  To  determine  necessary  and  sufficient  conditions  for  robust 
stabilization  and  (ii)  to  determine  whether  there  is  a  controller  that 
satisfies  these  conditions.  The  first  part  was  solved  by  Doyle  and  Stein 
(1981),  Chen  and  Desoer  (1982)  and  Vidyasagar  and  Klmura  (1984),  while  the 
second  part  was  solved  by  Vidyasagar  and  Kimura  (1984).  These  results  form 
the  point  of  departure  for  Task  3,  and  are  summarized  below. 

To  make  the  theorem  statements  more  compact,  we  introduce  three 
classes  of  uncertainty.  In  each  class,  r  is  a  given  function  belonging  to 
the  set  S  of  stable  transfer  functions.  Let  Pq  be  a  given  plant  descrip¬ 
tion,  not  necessarily  stable,  but  free  of  jw-axis  poles.  Then  the  class 
A( Pq , r )  of  additive  perturbations  consists  of  all  plants  P  that  have  the 
same  number  of  open  RHP  poles  as  Pq,  and  satisfy 

||  P  (jw)  -  Po(jw) | |<  |r  (jw)  |  for  all  w.  (3.1) 

The  class  M(Pq,t)  of  multiplicative  perturbations  consists  of  all  P  of  the 
form  (I+M)Pq,  such  that  P  has  the  same  number  of  open  RHP  poles  as  Pq,  and 


in  addition. 


i 


|  |M(jw) j |  <  |r  (jw) I  for  all  w 


(3.2) 


Finally,  set  (Ng  Dg  )  be  a  particular  r.c.f.  of  Pg.  Then  as  discussed  In 
Task  1,  the  class  S(Pg,r  )  of  stable  factor  perturbations  consists  of  all 
plants  P  of  the  form  N(s)  [D(s)J  ,  where 

IN  (jw)  -  N0  (jw)|| 

{  <  |  r(jw)  |  for  all  w.  (3.3) 

D  (jw)  -  Dg  (jw)« 

In  the  case  of  stable  factor  perturbations,  there  are  no  restrictions  on 
the  relative  numbers  of  RHP  poles  of  P  and  of  Pg 

Theorem  3.1  (Doyle  and  Stein  1981,  Chen  and  Desoer  1982) 

A  controller  C  that  stabilizes  Pg  also  stabilizes  all  plants  in  the  class 
A  (Pg,  r  )  if  and  only  if 

|  | [ C  (I  +  Pg  C)  -1  ]  (jw)  ||.  jr  (jw)  |  <  1  for  all  w.  (3.4) 

C  stabilizes  all  plants  In  the  class  M  (Pg,  r)  if  and  only  if 

| | [Pg  C  (I  +  Pg  C  )  -1  ]  (jw)  ||.  |r  (jw)  I  <  1  for  all  w.  (3.5) 

Theorem  3.2  (Vidyasagar  and  Kimura  1984).  Suppose  C  stabilizes  Pg^ 

and  choose  an  &.c.p.  (A,B)  of  C  such  that 

A  (s)  Dg  (s)  +  B  (s)  Ng  (s)  =  I,  for  all  S.  (3.6) 

Then  C  stabilizes  all  plants  in  the  class  S  (pg,  r)  if  and  only  if 

||([  A  B]  r)  (jw)  ||  <  1  for  all  w.  (3.7) 

Theorem  3.3  (Vidyasagar  and  Kimura  1984).  Let  (Ng,Dg),  (Dg,  Ng) 
be  a  particular  r.c.f.  and  a  particular  ft.c.f.  of  Pg,  and  choose 
X,Y,X,Y  in  M(S)  such  that 

X  Ng  +  Y  Dg  »  I,  Ng  X  +  Dg  Y  »  I  .  (3.8) 


38  - 


Then  there  exists  a  controller  C  that  stabilizes  all  plants  in  the  class 
A  (Pg,r)  if  and  only  if  there  exists  an  ReM  (S)  such  that 

||[(X"D+  DRD)  r]  (jw)  | |  <  1  for  all  w.  (3.9) 

There  exists  a  controller  C  that  stabilizes  all  plants  in  the  class 
M(Pg,r)  if  and  only  if  there  exists  an  ReM  (S)  such  that 

| | [ (N  X  +  NRD)  r  ]  (jw)  | |  <  1  for  all  w.  (3.10) 

There  exists  a  controller  C  that  stabilizes  all  plants  in  the  class 
S(P0,r)  if  and  only  if  there  exists  an  ReM  (S)  such  that 

| | r  (l  Y  X]  +  R  [-N0  D0  ])  (jw)  | |  <  1  for  all  w.  (3.11) 

Now  we  come  to  Task  3.  The  objective  of  this  task  is  to  determine 
conditions  for  the  existence  of  robust  decentralized  controllers,  in  the 
special  case  of  a  two-channel  system.  We  restrict  attention  to  the  case 
where  Pg  is  stable.  In  this  case  the  problem  is  one  of  determining  whether 
or  not  there  exists  a  decentralized  controller  C  that  satisfies  (3.4), 

(3.5)  or  (3.6),  as  appropriate.  It  turns  out  that  this  can  be  translated 
into  a  restriction  on  R  using  Theorem  2.8. 

Theorem  3.4  Suppose  Pg  is  2x2  and  decentrally  strongly  stabilizable. 
Let  Ci  be  a  diagonal  controller  that  stabilizies  Pg,  and  let  Pj  =  Pg 
(I  +  Cj  Pg  )  “1.  Let  qeS  be  given.  Then  there  exists  a  decentralized 
controller  that  stabilizes  all  plants  in  the  class  A  (Pg,q  )  if  and  only 
if  there  exists  a  matrix  R  s2x2  such  that 

r12  =  ~Cpl > 12  lRl»  r  21  “  _  (pl>  21  lRl>  (3.12) 

I  I qC !  (I-Pi  Cj  )  +  q  (I-  Pi  q  )  ||„  <  1.  (3.13) 
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There  exists  a  decentralized  controller  that  stabilizes  all  plants  in  the 


class  M(Pg,q)  if  and  only  if  there  exists  a  matrix  ReS2x2  satisfying  (3.12) 
such  that 

||q  Pt  CX  +  q  Pj  R(I  -  Pj  CL  )  ||  .  <  1.  (3.14) 

Finally,  let  (Nq.Dq)  be  an  r.c.f.  of  Pq  and  let  U  be  the  unique  unimodular 
matrix  in  s2x2  such  that 


N0 

-00- 


lJ  = 


Pi 

a 

-I-  Ci  Pi_ 


(3.15) 


Then  there  exists  a  decentralized  controller  that  stabilizes  all  plants  in 
the  class  S(PQ,r)  if  and  only  if  there  exists  a  matrix  R  s2x2  satisfying 
(3.12)  such  that 


||  (  U  [  I  Ct  ]  +  R  [  -  PL  I-  Pi  q  1)  q  ||.  <  1 


(3.16) 
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5.  TASK  4 


OPTIMAL  DECENTRALIZED  FILTERING  FOR  A  TWO-CHANNEL  SYSTEM 

In  this  task,  we  propose  an  iterative  solution  procedure  based  on 
Theorem  2.8.  Suppose  P  is  2x2,  W  is  a  given  2x2  weighting  matrix,  and  we 
want  to  minimize  the  cost  function 

J  -  | |  W  (I  +  PC)  “1  |  |„,  (4.1) 

where  the  norm  |  |  •  |  |oo  is  defined  by 


||f||o»  =  ess.  sup  I  I F  (jw)||,  (4.2) 

w 

and  we  want  to  find  the  optimal  decentralized  C.  If  P  is  decentrally 
strongly  stabilizable,  then  Theorem  2.4  allows  us  to  assume  without  loss  of 
generality  that  P  is  stable.  In  this  case,  the  set  of  all  (I+PC)“1  that 
result  from  a  stabilizinbg  controller  is  given  by  I-PR,  as  R  varies  over 
M(S).  However,  only  certain  choices  of  R  lead  to  decentralized  control¬ 
lers,  as  shown  in  Theorem  2.8.  Thus  the  problem  is  one  of  minimizing 
|  |w( I-PR)  |  |oo,  subject  to  the  constraint  that 

r12  =  "  P12  lRl>  r21  =  “  P21  lRl*  (4.3) 

Let  us  denote  R  by  y ;  then  (4.3)  implies  that 


r12  =  -  P12  Y»  r21  =  ~P21  Y  i  rn  r22  “  Y  (1"P12  P21  Y  )•  (4.4) 


The  cos t  function  to  be  minimized  is | |W( I-PR) | |»,  where 


I-PR  =■ 


1-  (Pll  rll  +  Pi 2  r2 1  ) 
r  (P21  rll  +  P22  r2 1  ) 


-  (Pll  rl2  +  P12  r22  )~ 
1-  (P21  r12  +  P22  r22  L 


-  Y  P12  P21  ~  Pi 1  rli  -  Pi i  Y  “  P12  r22 

Z  P21  rll  “  P22  P21  Y  I"  Y  P12  P21  “  P22  r22 

The  norm  of  this  quantity,  which  is  an  af f ine  function  of  the  vector 
[y  rji  r22  1*  Is  to  he  minimized  subject  to  the  constraint 

rll  r22  “  Y  (1-  P12  P21  Y  )  °  0.  (4.6) 

This  can  be  done  using  nonlinear  programming. 
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CONCLUSIONS 


With  various  degrees  of  detail  and  complexity  we  have  obtained  signifi¬ 
cant  results  on  all  of  the  proposed  tasks,  plus  an  additional  subtask 
together  with  two  robust  controller  designs  for  a  flexible  structure  model. 
More  specifically,  two  numerical  procedures  were  given  for  construction  of 
two  different  types  of  robust  controllers.  In  addition,  the  concept  of 
fixed  modes  was  extended  to  that  of  fixed  determinant,  which  is  important 
for  decentralized  control.  Furthermore,  it  was  demonstrated  that  for  a 
stable  system,  then  exists  a  class  of  robust  decentralized  controllers  for 
a  2x2  system.  Finally,  an  iterative  approach  was  proposed  to  solve  optimal 
decentralized  filtering  for  a  two  channel  system.  It  appears  evident  that 
in  this  Phase  I  preliminary  effort  we  have  demonstrated  the  feasibility  of 
the  approach. 
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APPENDIX  A 
(Proofs) 


In  this  appendix  we  give  the  proofs  of  all  results  that  are  original  to 
this  proposal. 

Proof  of  Theorem  2.2  First,  let  K  be  a  block-diagonal  matrix.  Then, 

( I ,K)  is  an  Jt.c.f.  of  K,  so  that  $(P,K)-|l  D4-K  N|=*|D+KN|.  Now 

|d+kn|»  |d|  |i+kp|  (a.i) 

Expanding  the  determinant  |l+KP|  using  the  well-known  formula  for  the 
determinant  of  the  sum  of  a  diagonal  matrix  and  an  arbitrary  matrix  (see 
e.g.  Vidyasagar  1985)  gives 

<p  (  P,  K  )  =  |d|.  1  +  E  principal  minors  of  KP  (A. 2) 

Now  the  minors  of  KP  can  be  expanded  by  the  Binet-Cauchy  formula  (Vidyasa¬ 
gar  1985 ) >  This  gives 

(KP)  I#I  -  E  Kt,J  PJ»I  <A*3> 

<(»  (P,K)  -  |D|  +  E  Kj,J  pJ»I  l°l 

However,  (Vidyasagar  1985)  the  quantity  Pj,i  |D|can  be  expressed 
as  a  minor  of  F,  and  several  minors  Kj,j  are  structurally  equal  to  zero 
because  of  the  block-diagonal  nature  of  K.  Thus  the  only  minors  that 
appear  in  (A. A)  are  the  minors  of  F  referred  to  in  the  statement  of  Theorem 
2.2.  If  A  denotes  the  g.c.d.  of  these  minors,  then  each  term  in  (A. 4)  is 
a  multiple  of  £,  whatever  be  K.  Hence  c  divides  <{>(P,K)  for  all  K,  and 
finally  divides  f,  which  is  the  g.c.d.  of  all  $(P,K). 


To  complete  the  proof,  we  show  that  f  divides  c.  First,  let  K*0. 

Then,  from  (A. 4),  <t>(P,K)=  |D  | .  Hence,  from  (2.3),  f  divides  |d|.  Next,  let 
K  be  a  matrix  with  a  single  nonzero  element,  which  equals  one.  Then  the 
summation  in  (A. 4)  reduces  to  a  single  minor  of  F,  wherein  one  row  of  D  is 

replaced  by  a  row  of  N.  Thus  f  divides  this  minor.  (This  part  of  the 

argument  is  very  similar  to  that  in  Vidyasagar  and  Viswanadham  (1983)  and 
is  therefore  only  sketched.)  Next  we  choose  K  to  contain  two  one's  and  the 
rest  zero.  This  will  show  that  f  divides  appropriate  minors  of  F  con¬ 
taining  two  rows  of  N  and  the  rest  from  D.  Finally,  f  divides  the  g.c.d. 
of  all  these  minors,  which  is  c, . 

Proof  of  Theorem  2.3  Clearly  g  divides  f  since  g  is  g.c.d.  of  a  larger 
set.  To  prove  that  f  divides  g,  write 

<j>  (P,C)  =  |l  +  PC  | .  |DC  |.  |D.  | 

*  |DC I .  |D  |  +  E  minors  of  P.  |d|.  minors  of  C.  |dc  | 

=  |DC  | .  | D |  +  E  minors  of  F.  minors  of  [  Dc  Nc  ]  (A. 5) 

However,  every  function  in  (A. 5)  is  a  multiple  of  C=f;  hence  <j>(P,C)  is  a 
multiple  of  f.  Finally,  the  g.c.d.  of  all  such  $(P,C)  is  a  multiple  of  f, 
i.e.  f  divides  g. 

Proof  of  Theorem  2.4  In  (Vidyasagar  1985,  Theorem  (5.3.10)),  it  is 
shown  that  if  C  stabilizes  P,  and  Cj  stabilizes  Pj,  then  C+Cj  stabilizes 
P.  Now,  if  C  and  Cj  are  both  decentralized,  so  is  C+Cj.  This  proves  that 
C+DS(Pj)  is  a  subset  of  DS(P).  Next,  suppose  C  is  decentralized  but 


unstable.  Then,  from  the  proof  of  the  above-cited  theorem,  we  know  that, 
for  sufficiently  small  a,  the  controller  aC  does  not  belong  to  S(Pj)  (and 
therefore  certainly  not  to  DS(Pj)),  even  though  (l+a)CeS(P).  Since  C  is 
decentralized,  so  is  (l+a)C.  Thus,  (l+a)C  DS(P)  but  otC  DS(Pi),  i.e. 
C+DS(Pi)  is  a  strict  subset  of  DS(P). 

Proof  of  Theorem  2.5  Suppose  C  is  stable  and  diagonal.  Then,  letting 
denote  the  ii-th  element  of  C,  we  get 

|D  +  CN|  =  |  dt  +  C!  nj  ...  dm  +  cm  nm  | 

=  |D  |  +  E  a  minor  of  c  .ua 
a  *  o 

»  | D |  +  a  multiple  of  g.c.d.  ua  (A. 6) 

a  *  o 

In  order  for  this  quantity  to  be  a  unit,  |d|  and  g.c.d.  va  mu8t  satisfy  the 
parity  interlacing  property. 

Proof  of  Lemma  2.6  see  Vidyasagar  et  al.  (1985). 

Proof  of  Theorem  2.8  Let  C-Cl-RP)"^.  Then  C=Y_1X,  where  Y=*I-RP,  X=R. 
Let  the  subscript  i  denote  the  i-th  column  of  a  matrix.  Then,  by  Cramer's 
rule , 

c12  “  Oi  c\2  *  o  if  and  only  if  yj  A  xi  “  0  ,y 2  A  X2  “  0  •  (A. 7) 

Let  ei  denote  the  i-th  elementary  unit  vector,  i.e.  the  i-th  column  of  the 
identity  matrix.  Then  y^  *ei~Rpi,  y2“e2-®P2*  Similarly,  x^*rj,  X2=*f2* 


yj  A  xj  =  0  if  and  only  if  (cj  -  RPj  )  A  rj  =0 


(A. 8 


But  Rpl-rj  pH  +  r2  P2l,  Hence 


yi  A  xj  =0  if  and  only  if  cj  A  rj  -  (rj  Pi i  +  r2  p2i  )A  rj  =  0 

if  and  only  if  cj  A  ti  -  p2l  r2  A  ri  =  0  (A. 9 

since  rj  A  =0.  Now 


el  A  rl 


1  rn 

0  r21 


r2 1 «  r2  A  rj  =  rj  A  r2  =  -|r|. 


(A. 


Hence  (A. 9)  becomes 


y2  A  xi  =  0  if  and  only  if  r2i  =  -  p2i  |r|.  (A.l 

Similarly, 

y2  A  X2  =  0  if  and  only  if  r[2  =  -P12  lRl 
Proof  of  Theorem  3.4  As  shown  in  (Vidyasagar  1985,  p. 
choose 
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as  a  doubly  coprime  factorization  of  Pq.  Moreover,  C+Ci  stabilizes  Pq  if 
and  only  if 


(A.l 

) ,  one  can 


Ci  =  R  (  I  -  P1R  )  -1  t  R  e  s2x2  . 


(A. 
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0  OOOOD+OO 

0.  OOOOD+OO 

o.  ooo'.: 

-9  0139D-01 

3.  0395D+02 

0.  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOOD+OO 

0.  000 CD 

2  0721D+04 

-1. 1 988D+07 

0  OOOOD+OO 

0  OOOOD+OO 

0.  OOOOD+OO 

0  0000" 

-7  7263D+03 

1.  81  50D+06 

0  OOOOD+OO 

0. OOOOD+OO 

0.  OOOOD+OO 

c.  ooc\ : 

3  5874D+02 

-1.  0678D+05 

0  OOOOD+OO 

0  OOOOD+OO 

0.  OOOOD+OO 

o.  oooo : 

1B842D+02 

-3.  7800D+04 

0.  OOOOD+OO 

0.  OOOOD  +00 

0.  OOOOD+OO 

0.  000  jL' 

-5  2522D+01 

1  21 01D+04 

0  OOOOD+OO 

0.  OOOOD+OO 

0  OOOOD+OO 

0.  OOOOL 

-2. 0237D+00 

2.  2510D+01 

0.  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOO'L 

1  5532D+00 

-1. B653D+02 

0  OOOOD+OO 

0  OOOOD+OO 

0  OOOOD+OO 

0  OOOOD 

-2  1345D+00 

5.  9849D+01 

-8.  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOOD 

0  OOOOD+OO 

0  OOOOD+OO 

1  OOOOD+OO 

-8. OOOOD+OO 

0  OOOOD+OO 

0  000  Di 

0  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOOD+OO 

1. OOOOD+OO 

-B.  OOOOD  +  OO 

0.  OOOOL 

0  OOOOD+OO 

0  OOOOD+OO 

0  OOOOD+OO 

0  OOOOD+OO 

1 .  OOOOD+OO 

-8  OOOOL 

0  OOOOD+OO 

0.  OOOOD+OO 

0.  OOOOD+OO 

0. OOOOD+OO 

0.  OOOOD+OO 

1 .  OOOOL 

Columns  19  to  19 
- 1  0486D+07 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0.  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
0  OOOOD+OO 
-8  OOOOD+OO 


-  52 


4. 6896D+03  4.  02350+05 

3  51080+02  -0.  2540D+04 

1  04540+01  2.  3600D+04 

9  1677D-01  1 . 28 1 80+02 

1  8140D-02  -7.  76730+01 

2  1341 D— 03  -1.  3664P+00 

0  0000D+00  -4.  50690-01 

1  •7169D+01  1 . *9560+04 

2.  5557D+00  -2  5S50D+03 

4,  5383D  -02  1  56420+02 

fcj  3230D-02  5.  25890+01 

-i  6140D-02  -1.71900+01 

-2  1341D-03  5.  534)60-02 

0  G000D+00  2.  76590-01 

-2. 1341D-03  0. OOOOD+OO 

0  0000D+00  0.  00000+00 

0. OOOOD+OO  0.  OOOOD+OO 
O  OOOOD+OO  O.  OOOOD+OO 

0  OOOOD+OO  O.  OOOOD+OO 


Columns  1  to  6 

3  5000D+01  -6.  9516D+02  1.59380+04  -31472D+05  53675D+06 

Columns  7  to  12 

2  6991D+08  3.85000+01  -1 .  07450+03  1.9195D+04  -6.  2458D+05 

Columns  13  to  18 

-6  7300D+07  3.  0853D+08  8.  00000+01  -2  56000+04  4  0960D+05 

Columns  19  to  19 
1  0486D+07 


-6  9M7D+03  0  OOOOD+OO 


-2  8090D  *-07 


6  31560+06 


-3  2768D+06 


-  53  - 
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